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o\ , 

, We present a new method to study the properties of heavy quarks at finite 

temperature. It combines non-relativistic QCD with an improved gluonic 
action on anisotropic lattices. The efficiency of the approach is demon- 
strated by the first non-perturbative calculation of the temperature de- 
If^ I pendence of low- lying quarkonium "pole" masses. For ground state meson 

masses in the region between charmonium and bottomonium we find only 
very little variation up to our highest temperature which corresponds to 
^ ' T w 1.2 Tc while first excited states indicate a large mass shift. 

(N ■ 

O 

l> . 

O ■ 1 Introduction 

l> 

^ ' At high temperature hadronic matter is expected to undergo a phase transition 



to the quark-gluon plasma (QGP). In recent years there have been significant 



D . theoretical and experimental developments to study it's properties. High en- 

^ ' ergy heavy-ion experiments designed to detect spectral changes of hadrons in 

hot media have already been started. In the search for possible signals of the 
^ ' QGP heavy quarkonium states are among the simplest probes that allow to 

■ test the structure of the QCD vacuum. For instance J/\l/-suppression due to 

colour screening has been proposed to probe deconfinement Quarkonium 
production is rather well understood in hadron-hadron collisions 0. How- 
ever, the behaviour of bound states of heavy quarks in a strongly interacting 
medium close to the deconfinement temperature Tc is still largely uncertain. 

Previous theoretical results on the temperature dependence of the masses of 
the rjc, the J/\l/, and the are not yet completely satisfactory. Different calcu- 
lations give model dependent results. Some predict a significant decrease [Q-]^ 
while others suggest that the masses stay constant M or even rise with temper- 



ature PHT0|. The assumptions for instance about a temperature dependence 
of the string tension, cr(T), the strong coupling constant, dsiT), the effective 
quark mass Mq{T) and the gluon condensate {Q\F^^F'^''\Q) = G'^{T) primar- 
ily depend on perturbation theory and may not be reliable in the temperature 
region under consideration. 
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Progress in numerical simulations of QCD makes reliable predictions about 
hadronic properties in the non-perturbative regime possible. Recently a non- 
relativistic approximation of QCD has been used to reproduce the experimen- 
tal spectrum for heavy quarkonia with high precision ||TI| , p!2| . 



At finite temperature a complication arises because Lorentz invariance is ex- 
plicitly broken so that Green's functions defined by correlators in Euclidean 
time and space directions are controlled by different phenomena. The heavy 
quark potential which is confining in the low temperature region becomes 
Debye-screened at high temperature. On the other hand the pseudo-potential 
from spatial Wilson loops is confining for all temperatures irsHTSl . So far most 
of the work has been concentrated on the calculation of screening masses ob- 
tained from spatial correlators which do not have a direct connection to the 
physical mass of a resonance as defined for instance by the position of a peak 
in the spectral function. In recent investigations of the temperature dependent 
structure in the light mesonic channels screening masses from spatial correla- 
tions have been found to differ from effective pole masses from correlations in 
the Euclidean time direction P,|To| . The mass shift of hadrons made of light 
quarks is expected to be mainly controlled by chiral symmetry restoration. 
A realistic simulation would require the inclusion of light dynamical fermions 
which is expensive. To begin with, it is advantageous to study the temperature 
dependence of the spectrum with heavy quarks. In this case the influence of 
light quarks is less important so that a quenched simulation is a reasonable ap- 
proximation. The binding energy for cc (~ 0.63 GeV) and bb (~ 1.1 GeV) 
is large compared to the deconfinement temperature ~ 150 — 250 MeV. At 
intermediate temperatures the spectral width of the low lying qq bound states 
is expected to be small so that there will still be a clear distinction between 
the continuum and the lowest resonance. The binding energy of quarkonia 
decreases with temperature. The dissociation temperature of the (T, T') is 
expected to be in the region of T ^ (2.6, 1.1) ^7^. However, these values 
still contain a model dependent uncertainty. 

Our approach is based on an improved gluonic action for the light degrees of 
freedom and a non-relativistic formalism (NRQCD) for the heavy quarks |18| ]. 
Finite temperature NRQCD (FT-NRQCD) uses anisotropic lattices [pJl-pl 



to achieve a finer resolution in time direction. A large number of Matsubara 
frequencies is necessary to accurately measure temporal meson propagators. 
The problem is first approached in quenched QCD by considering a quark- 
antiquark pair propagating in Euclidean time direction in a gluonic medium. 
To show the feasibility of the new approach and the significance of its results 
we calculate temporal meson correlators for bare quark masses between 2 and 
6 GeV. 
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2 Heavy meson spectrum at finite temperature 



Heavy quarkonia are small and tightly bound. We know that asymptotically 
for infinite quark mass a potential picture will give the correct description. 
The Cornell potential 

V{r) = ar-— (1) 



reproduces the experimental spectrum for charmonium and bottomonium 
quite well. In a thermal medium of temperature T > the potential is modified 
by colour screening which can be parameterized in the form 



cr Of 
Vrir) = -7^(1 - exp(-/i(T)r)) exp{-fi{T)r) 



(2) 



Eq. is equivalent to the Cornell form with a temperature dependent string 
tension 0"(T) = cr/{fi{T)r) (1 — exp(— /i(T)r)) and a screened Coulomb term 
as{T) = exp(— /i(T)r). From fig. ^ we see that Vrir) changes differently 
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Fig. 1. The heavy quark potential (a) as a function of the radius for 3 different 
temperatures and the terms contributing to the total energy together with the 
dissociation energy in GeV as a function of the Debye-mass for cc (b) and bb (c). 

with temperature for small r < 0.3 fm and large r > 0.4 fm. The temperature 
dependence is entirely contained in the Debye-mass. It has been argued 
that due to string breaking /^(T) will be different from zero even for T = 
and that the effective string tension does not vanish immediately above Tc 
modeling non-perturbative interactions in the plasma. Although the precise 
functional dependence is not very well known the string tension and the strong 
coupling will decrease when the temperature and thereby /i(T) is increased. 

First qualitative insight in the behaviour of QQ-sta.tes can be obtained from a 



semi-classical picture |23[. As a consequence of the uncertainty relation which 
forces {p'^){r'^) ~ 1 the kinetic energy, E^in = p'^/m ^ 1/mr^, decreases with 
the quark separation. Minimizing the total Energy, E{r) = 2m + Ekin + V, the 
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ground state is found where the decrease of the kinetic energy is compensated 
by an increase of the potential energy. The result is that the binding radius 
ro(/x) increases with n{T). A minimum of E{r) exists as long as the screening is 
not too strong. The meson will dissociate once the Debye-mass n{T) becomes 
larger than a critical value fic- Below this value the meson mass depends on 
the temperature. In general the sign and the magnitude of the mass shift, 
AM = M{T) — M{T = 0), depend on the details of the balance of potential 
(Coulomb and string) and kinetic energy as is shown in figs. |l|b and c. In the 
potential of eq. ^ masses can decrease or increase with T depending on the 
size of the QQ bound state. 

Understanding the basic mechanism in the semi-classical approximation we 
can move towards a quantitative understanding and calculate quarkonium 
wave-functions in a non-relativistic potential model. The mean quark veloc- 
ity can be computed from averages obtained by a numerical solution of the 
Schrodinger equation with a potential of the form given in eq. 0, 

, EiT)-aiT){r)T + a,iT){l/r)T 




Fig. 2. Dependence of the velocity, radius, mass and dissociation energy in the (IS) 
state on the Debye-mass for bb (a) and cc (b). The velocity and the meson mass 
have been rescaled to fit in the same frame. 

From fig. ^ we see that the average quark velocity decreases with temperature 
so that a non-relativistic ansatz for the heavy degrees of freedom seems to be 
justified for T > 0. 

However, beyond this simple potential picture we know that lattice data do 
not agree very well with the perturbative expectation for the form of the heavy 
quark potential at least up to T < 4 p4|-[26[]. A precise measurement of the 
Debye-mass is plagued by ambiguities in the form of the function used to fit 
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the heavy quark potential. 



Lattice simulations can provide a reliable determination of the meson spectrum 
without assumptions about the heavy quark potential. Thus, it is important 
to compare phenomenological and numerical results to gain further insight in 
the structure of the QCD vacuum at finite temperature. 

Compared to zero temperature additional care is necessary because any ex- 
citation acquires a finite lifetime at non-zero temperature P^. The spectral 



function p(j9, u) of an excitation with momentum p and energy u which is a S- 
function at T = will broaden. This effect will modify the meson propagator, 

Gm(t). 



3 FT-NRQCD 



We can go beyond a potential model approximation and calculate the action of 
heavy quarks in a gluonic heat-bath from first principles. The proposed method 
is based on the formalism of non-relativistic QCD [|l8l and a discretization of 



quantum field theories at finite temperature on anisotropic lattices [IS -21 



NRQCD is an effective theory that has been successfully applied to bottomo- 



nium and charmonium at zero temperature |TT,T^. It allows an efficient and 
accurate calculation of heavy quark propagators. At zero temperature the 
binding energy of a 66-system can be estimated from it's distance to the open 
b threshold, Eb = M^g— ~ 1.1 GeV. The binding energy is large compared 
to the deconfinement temperature, ~ 150 — 250 MeV, so that we expect a 
smooth evolution of the system away from it's ground state at T = as the 
temperature is increased. From previous considerations of potential models we 
expect the average quark velocity to decrease with temperature so that there 
will be a range of temperatures T > where FT-NRQCD is applicable. 

The total action for a system of heavy quarks in a gluonic heat-bath naturally 
splits into a relativistic part for the gluons and a non-relativistic term for the 
heavy quarks. 

S = Sr + Snr (4) 



The relativistic simulations use a tree-level improved gauge action [0,^] with 
a plaquette and a rectangle term on asymmetric lattices with separate spatial 
coupling jSa- and temporal coupling 
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The asymmetry parameter 7 = s^fi^jfi^ and the couphng /5 = ^f^T^ are 
defined in terms of the spatial couphng fifj and the temporal coupling fir- In 
principle, tadpole improvement can be implemented by factors 

= Re Tr a,)i/^ and Ur = Re Tr a,)i/^ 

which can be calculated from the average spatial and temporal plaquette as 
indicated or the mean link in Landau gauge. The anisotropy C, = aa/O'T > 1 
becomes equal to the asymmetry parameter 7 only in the continuum limit, 
/5 — > CX3. At finite gauge coupling, they differ by a renormalization factor t] = 
^/'j which can be determined nonperturbatively in a calibration procedure. 
The additional parameter 7 allows to have a large number of lattice points in 
time direction while keeping the temperature T = (A^,- a,-)^^ = {N^- ^ Oo-)^^ 
fixed. 

The non-relativistic action is derived from the Dirac equation by a Foldy- 
Wouthuysen transformation. The resulting effective field theory called NRQCD 
approximates relativistic QCD at small energies. Relativistic heavy-quark mo- 
menta are excluded from the theory by choosing a spatial lattice spacing 
~ Mq. The 1/Mq expansion underlying the formalism of NRQCD will 
hold as long as the average quark velocity is small, f <^ 1. The non-relativistic 
action has the form 

SNR = ^HDt + Ho + 6H)^ (6) 

where Dt denotes the covariant time derivative, Hq is the kinetic energy op- 
erator and 6H is the leading relativistic and finite-lattice-spacing correction. 
Here we include all spin-independent relativistic corrections to order Mqv^ and 
spin-dependent corrections to order Mqv^. Modifications to the corresponding 
form of Snr for T = as given in ref. [^] appear only in the correction 6H 
where the improved temporal derivative and the chromoelectric field strength 
introduce additional factors of ^. To simplify matters, tadpole improvement 
of Snr was implemented only for spatial links. In our case the spatial lattice 
spacing is considerably larger than the temporal spacing so that the mean 
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temporal link Ur is very close to unity ^1 . 

The heavy quark propagators are computed using the evolution equation 



Gt+a-r 



° ^ f/i ( 1 



arSH) Gt 



(7) 



The non-relativistic quark propagator is not periodic in time and can be eval- 
uated at times larger than Nt-/2. A symmetric (antisymmetric) propagator 
can be constructed by explicitly adding (subtracting) the contributions from 
mirror charges. On a Euclidean lattice at temperature T = 1/(0,- A^^) ther- 
mal Green's functions can be evaluated only on a discrete set of frequencies 
ojn = 27inT, n = l...A^^. It is obvious that a good resolution requires a large 
number of grid points in temperature direction. Thermal meson states 



GUx 



)r(-'=)(r) J2Gt{x-s)T^''\s) 



corresponding to ^5*1 and ^S'o-states were constructed from the quark propa- 
gators. Interpolating operators of the form T{r) = fi*^^*" $(r) with 11 different 
combinations of local and smeared trial wave-functions at the source and sink 
were used. The wave-functions were determined solving the Schrodinger equa- 



tion for a Breit-Fermi potential BT| for each value of the bare quark mass. 



The same set of wave-functions was used for all temperatures. 



4 Simulation details and results 



For the relativistic simulations we choose an asymmetry parameter 7 = 4 so 
that our largest lattice size, 16^ x 64 corresponds to a symmetric lattice at zero 
temperature. The temperature was varied by reducing the value of from 
64 to 24 and 16 for fixed spatial lattice size, N^r = 16. The value of the bare 
coupling, P = 4.31466, corresponds to the critical coupling Pwnson = 5.8941 
P^ , |S3[] for 7 = 1 on a lattice with A^^nticai = 5_ this way it is possible 
to get a rough estimate of the spatial lattice spacing from the string tension 
^/o^a = 0.2734(37) at Pwnson = 5.8941 ||33|. We expect the inverse spatial 
lattice spacing to be in the region = ^/a/0.2734: 1.5 GeV. However, a 
value of 7 7^ 1 will modify this correspondence and change the value of the 
spatial lattice spacing. The remaining free parameter, the bare quark mass, 
was varied in the range a^Mq = 1.5,2,2.5,3,3.5,4. After generating ^ 100 
independent gauge field configurations we determined zero momentum meson 
propagators for ^Sq and ^Si states. For local ^Sq states we also measured finite 
momentum propagators. 
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Before we can measure the meson spectrum we have to determine the asym- 
metry ^ = ttrj/ar and set the scale a^. The heavy quark potential at T = can 
be measured either from spatial or temporal Wilson loops. The asymmetry 
parameter is determined in a calibration procedure from spatial and temporal 
potential differences. The parameter ^ can be calculated most accurately at 
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Fig. 3. The ratio of spatial and temporal potential differences (a) and heavy quark 
potential from temporal and spatial Wilson loops (b). 



small distances where the statistical error is small. In fig. ^ we see that for 
the smallest separation R = 2 the ratio of potential differences agrees with 
^ = 4.65 for all separations T > 2. All values for larger separations R = 2.83 
and -R = 3 are in accord with this value within the statistical errors. Fig. 
shows that spatial and temporal Wilson loops give the same physical potential 
when this value of the asymmetry parameter is used. Knowing ^ and a^^y/a 
an estimate of the temperature can be obtained from the relation 



T 



T 

7^ z 



19 

— ^0.8,1.2 (iV, = 24,16) 



where we use a value of Td \fo = 0.625 A least squares fit to the static 
potential with the functional form V{R) = Vq — it / {12R) + a R gives the result 
a^cT = 0.16. This translates to ^ 1.1 GeV when the phenomenological 
value for the string tension, y/a = 427 MeV, is used. To show the consistency a 
second estimate for the spatial lattice spacing is determined from the hyperfine 
splitting of the ground state, AE^yp. Two conditions are needed to fix the 
unknown parameters and Mq. First we set the scale and determine 
from Mkin = Mjj^ = Mx — AE^yp where we use the experimental value for 
the mass of the T{1S) = 9.46037(21) GeV The kinetic mass Mk^n is 
determined from a fit to a non-relativistic dispersion relation, 
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Table 1 

The kinetic mass and hyperfine splitting of the ground state for various values of 
the bare quark mass for ^ = 4.65 and AEhyp between 15 and 25 MeV. 







a-J = 




dr^Ehyp 






~ e a-^ 


1.5 


15.4 ( 2) 


2.85 


GeV 


0.0011 - 


0.0019 


2.0 


20.1 ( 3) 


2.18 


GeV 


0.0015 - 


0.0025 


2.5 


24.8 ( 4) 


1.77 


GeV 


0.0018 - 


0.0030 


3.0 


29.6 ( 5) 


1.48 


GeV 


0.0022 - 


0.0036 


3.5 


34.5 ( 7) 


1.27 


GeV 


0.0025 - 


0.0042 


4.0 


39.8 (10) 


1.10 


GeV 


0.0029 - 


0.0049 



where M{p) is the mass extracted from finite momentum propagators. From 
previous calculations AE^yp is expected to be in the region between 15 and 
25 MeV at our presumably large value of between 1.1 and 1.5 GeV [Bl . 
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Fig. 4. The kinetic mass of the ^5*0 meson as a function of the bare quark mass. 



The value of = ^ M^^/(,^Mfcj„) is used to compare the expected AE^yp with 
the measured ^5*1 — ^Sq splitting. From figs. | we see that the best agreement 
is achieved for Mq = 3.5 which corresponds to ~ 1.3 GeV. There is a 15% 
difference between this value and the estimate from the string tension. Both 
values are consistent within the expected accuracy. A discrepancy of this size 
has also been found in NRQCD studies at T = A more stringent test 
and a better determination of a„ will be possible once the 15" — IP splitting 
has been measured . 
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The easiest quantities to calculate are ^5*0 propagators for point sources. Scaled 
propagators, H(T) = G(T)/G(T = 0), show directly the changes in the spec- 
trum due to the temperature. The zero point energy which in principal can be 
computed using weak coupling perturbation theory is the same for all three 
temperatures because we vary only Nr keeping all other parameters fixed. 
The mass shift AM(T) = MmesoniT) — MmesoniT = 0) can be extracted from 
the slope of the scaled propagator at large time steps. A positive slope indi- 
cates an increase while a negative slope will be seen if masses decrease with 
temperature. 

Our numerical results show that the signal for QQ bound states persists to 
high temperature, T ^ 1.2 T^. Compared to T = the ground state meson 
propagators change very little. The observed temperature dependence of the 
scaled propagator (fig. ^ is weak. Looking at it in more detail we see that the 
ground state meson mass decreases with T. Fig. | shows the scaled propagators 
for Nr = and 24. The scale is logarithmic and we expect a linear form for an 
exponential behaviour, G{t) oc exp (—AM t). For A^,- = 24 which corresponds 
to T ^ 0.8 Tc the scaled propagator agrees with unity. It does not show any 
significant change with temperature in the entire mass range, 1.5 < a^Mq < 4. 
For A^^ = 16 which corresponds to a temperature T ^ 1.2 Tc we observe a 
clear increase of the scaled propagator with Euclidean time. This signals a 
broadening of the spectral function and is an indication of a possible decrease 
of the effective meson mass at this value of the temperature. As expected, 
the effect becomes weaker as we go to larger quark mass. To illustrate the 
order of magnitude we included a straight line in fig. ^ which corresponds to a 
small shift of —AM = a~^(ln 1.06)/6 ^ 12 MeV. It is important to note that 
we used point sources and sinks which do not distort the spectral function. 
The observed effect is a signal that the mixture of states excited by the local 
source gets broader and possibly lighter with increasing temperature. If there 
is no cancellation in the sense that the ground state gets heavier while higher 
excitation become lighter then particularly the ground state mass will decrease 
with temperature. 



It is expected that first excited states dissociate immediately above the phase 
transition. In fig. ^ we see no effect at T 0.8 but a strong increase of 
the propagators at T 1.2 Tc. Compared to the ground state the increase is 
much stronger so that it seems hkely that the first excited state is dissociated 
at this temperature. Again, the order of magnitude of the effect is illustrated 
in fig. 1^ by a straight line which now corresponds to a shift of —AM = 
a~^(ln3)/10 ~ 240 MeV. Since there is no local operator projecting on the 
first excited state this result has been obtained with a smeared source and 
sink. The trial wave-function can only approximate the temperature wave- 
function causing a possible residual distortion of the quarkonium propagator. 
To overcome this problem a measurement of the temperature wave function 
is desirable [BB . 
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Our numerical results can be contrasted with an expectation from the den- 
sity matrix formalism. If we assume that the thermal meson state consists of 
a mixture of zero temperature states with a statistical weight given by the 
Boltzmann factor G(t,T) = exp (— Mj/T) exp (— Mj t) then we would ex- 
pect a decrease of the scaled propagator with t because more and more heavier 
states contribute as the temperature is increased. The observed increase of the 
scaled propagators cannot be explained by a simple temperature dependent 
admixture of higher mesonic states with negligible spectral width. 



5 Summary 



We have performed a first non-perturbative study of the heavy meson spec- 
trum and dispersion relation at finite temperature using FT-NRQCD. We 
found no change in the spectrum for temperatures below the deconfinement 
transition and strong non-trivial effects above T^. 

A comparison of S-state meson propagators at T = and T ^ 0.8 showed 
no significant differences in the entire Euclidean time range, t < 16. The local 
operator used for the ground state does not distort or re-weight the spectral 
function. At low time steps it gets contributions not only from the ground 
state but also from higher excitations. If this quantity is unchanged at large t 
so is the true ground state. 

A further increase of the temperature to T ^ 1-2 Tc showed a clear signal for a 
broadening of the spectral function of the ground state. This can be interpreted 
as an indication for a small decrease of the meson mass with temperature. 
The phenomenon is stronger the lighter the quark mass is. As expected we see 
that the first excited S-state is more susceptible to the temperature. Across 
the deconfinement transition we observe a strong change between T ^ 0.8 Tc 
and T ^ 1.2 Tc. This can be taken as a first indication for the dissociation of 
these states. 

Qualitatively all our observations are in accord with the expectations from the 
phenomenological picture of Debye-screening in a potential model. The heavier 
the mesons the smaller they are and the less they feel the screening. This has 
important implications for the deconfinement transition. The confined degrees 
of freedom are not released at a single temperature Tc but the transition to a 
quark gluon plasma proceeds in steps. The determination of the dissociation 
temperature of different bound states of heavy quarks is an interesting problem 



which deserves further efforts 36 



There is a discrepancy with other models that predict a strong effect already 
below T < Tc such as a weakening of cr(T). The string tension surely has 
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to vanish at T = Tc- However, this effect is a long range phenomenon. At 
smaller distances remaining strong interactions can still give rise to an effective 
string tension a{R,T) > 0. Furthermore, we see only a small change of the 
ground state mass up to T = 1.2 unlike predicted by most phenomenological 
models. A crossover between a behaviour dominated by the strong coupling 
constant at small distances and the string tension at larger separations is 
not observed in the accessible mass range. This gives further support to the 
observation that the perturbative form of the heavy quark potential is wrong 
in the temperature range under consideration. 

A precise measurement of the shift of the meson mass at higher temperature 



requires further work |^6[. For the future we intend to measure the wave- 
function and use it to improve the ground-state overlap of our operators. We 
expect that a multi-state fit with the improved operators allows a determina- 
tion of the spectrum at higher temperature. 
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Fig. 5. Local masses In (Gm(t)/G,„(t + a,-)) for the — ^Si splitting in units of 
from local-local (11), smeared- local (si) and smeared-smeared (ss) propagators 
compared with the expectation obtained from AE = 15 — 25 MeV. We start with 
the lowest bare quark mass a^Mq = 1.5 in the upper left corner and go in intervals 
of 0.5 to the lower right corner with a^Mq =4.0. 
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Fig. 6. Scaled local n=l ^Si meson propagators H{T) for 6 values of the bare quark 
mass starting with the lightest mass in the upper left corner on a logarithmic scale. 
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Fig. 7. Scaled smeared-smeared n=2 ^Si meson propagators H{T) for 6 values of the 
bare quark mass increasing from upper left to lower right corner on a logarithmic 
scale. 
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